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■ For the past 40 years, Brueckner theory has proven to be a most powerful tool 

to investigate systematically models for nuclear matter. I will give an overview 
of the work done on nuclear matter theory, starting with the simplest model and 
^ ' proceeding step by step to more sophisticated models by extending the degrees of 

Q , freedom and including relativity. The final results of a comprehensive hadronic 

. theory of nuclear matter are compared to the predictions by currently fashionable 

' two-nucleon force models. It turns out that a two-nucleon force can, indeed, re- 

' produce those results if the potential is nonlocal, since nonlocality is an inherent 

quality of the more fundamental fieldtheoretic approach. This nonlocality is crucial 
for creating sufficient nuclear binding. 



0^ 



1 Introduction 



> 

' By definition, nuclear matter refers to an infinite uniform system of nucleons inter- 

acting via the strong force without electromagnetic interactions. Tliis liypotlietical 
system is supposed to approximate conditions in tlie interior of a heavy nucleus. 
We shall assume equal neutron and proton densities; that is, we will consider sym- 
metric nuclear matter. This many-body system is charactarized by its energy per 
nucleon as a function of the particle density. Based upon various semi-empirical 
sources, nuclear matter is determined to saturate at a density po — 0.17 ± 0.02 
O \ iin~^ (equivalent to a Fermi momentum kp — 1.35 ± 0.05 fm"-'^) and an energy per 

2 ' nucleon f /A = -16 ± 1 MeV. 

Historically, the first nuclear niatter calculations were performed by Heisen- 
berg's student Hans Euler, in 1937.EJ This was just two years after Weizackercl had 
k>( \ suggested the semiempirical mass formula. Euler applied an attractive potential of 

5_j ■ Gaussian shape in second-order perturbation theory. 

\ Modern studies began in the early 1950's after a repulsive core in the nuclear 

potential had been conjecturedE It was obvious that conventional perturbation the- 
ory was inadequate to handle such singular potentials. Therefore, special methods 
had to lie_developed. This program was initiated by Keith Brueckner and co- 
workerguijO who applied — for the ground-state problem of nuclei — methods similar 
to those developed by WatsonEl for multiple scattering. Later, a formal basis for 
this new approach was provided by Goldstonecl who, using perturbation theoretical 
methods, established the so-called linked cluster expansion. The success of Brueck- 
ner theory in practical calculations stems from the fact that certain classes of linked 
diagrams can be summed in closed form up to infinite orders defining the so-called 



* INVITED TALK PRESENTED AT THE TENTH INTERNATIONAL CONFERENCE ON 
RECENT PROGRESS IN MANY-BODY THEORIES, SEPTEMBER 10-15, 1999, SEATTLE, 
WASHINGTON, USA; TO BE PUBLISHED IN ADVANCES IN QUANTUM MANY-BODY 
THEORY, VOL. 3 (WORLD SCIENTIFIC, SINGAPORE). 

t DEDICATED TO KEITH BRUECKNER ON THE OCCASION OF HIS 75TH BIRTHDAY. 



mbx: submitted to World Scientific on February 9, 2008 








-20 - 



I , I I I I , I , I I 

1 1.5 2 

kj, (fm"') 

Figure 1. Energy per nucleon in nuclear matter, £/A, versus the Fermi momentum, kp, applying 
the Argonne V14 potential. Solid curves represent Brueckner results, using the gap choice for 
the single particle potential, with the number of hole lines indicated in parentheses. The dashed 
curve is the Brueckner two hole-line result using the continuous choice. The two dotted lines are 
obtained in variational calculations. The box describes the area in which nuclear saturation is 
expected to occur empirically. (Dashed curve from Ref. all other curves from Ref. ^^.) 



reaction matrix G. All quantities are then formulated in terms of this G which — in 
contrast to the bare nuclear potential — is smooth and well behaved. 

The first numerical calculations ^applying Brueckner theory were perfcpmed 
in 1958 by Brueckner and Gammeu using the Gammel-Thaler potential£3 In 
the 19fiDaTn|-EJp4S Bethe and his students elaborated thpiiQiighly on Brueckner 
theory.ElrEjIlj'Efl Around 1970, elegant numerical methodaij'Ej for the soluticp^rLE 
the Brueckner equation were established and sytematic calculations performed. EZfE3 

In the mid 1970s, the nuclear many-body community was shaken by an ap- 
parent discrepancy between results from Brueckner Iheory and the variational ap- 
proch (the famous/infamous 'nuclear matter crisis'E3). This indicated that both 
many-body theories had to be reaxamined and more consistent calculations had 
to be perforttied. For Brueckner theory, the necessary work was done maialy by 
Ben Day.E3'Eil The variaticpal approach was |pursued by the Urbana group J23 the 
Rome-Pisa collaboration,c3 and many others.cJ As a result of this enormous work, 
quantitatively very close predictions were finally obtaped. from the two different 
many-body approaches using realistic NN potentials.c3E3E3 Based upon these com- 
prehensive investigations, it is now commonly accepted that both many-body ap- 
proaches are reliable for densities typical for conventional nuclear physics. 

Brueckner theory does not converge in powers of G, but in terms of the hole- 
line expapdon.Ej This is demonstrated in Fig. 1 for the case of the Argonne V14 
potential.E3 The solid lines in this figure are obtained from Brueckner theory using 
the conventional choice ('gap choice')c3 for the single particle potential; the num- 
ber of hole lines taken into account areJndicated in parentheses. The Brueckner 
results shown in Fig. 1 are from Dayrj they were recently confirmed by Baldo 
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Figure 2. Nuclear matter saturation as predicted by a variety of A'^A'^ potentials. Solid dots are 
the saturation points (i. e., minima of £/A versus kp curves) obtained in the Brueckner two hole- 
line approximation applying the gap choice for the single particle potential. Solid stars denote 
corresponding predictions that either include four hole-lines using the gap choice or two hole-lines 
using the continuous choice. The open star is the result by Brueckner and Gammel of 1958.® 



and coworkers.Ey The dotted lines in. Fig. 1 are the predictions by two variational 
calculations performed by Wiringa.cJ Within the uncertainties there appears to be 
agreement between Brueckner theory and the variational approach; particularly, 
around the minimum of the curves. Finally, the dashed curve in Fig. 1 is the lowest 
order Brueckner (i.e., two hole-line) result when the continuous choiceE3 is used 
for the single particle potentialJEl It is clearly seen that this latter result agrees 
well with the four hole-line Brueckner prediction using the gap choice. This may 
suggest that tVip rnritinons choicc is a more efficient way to perform Brueckner 
calculations.BEi3@ 

Systematic Brueckner calculatipns conducted with a variety of NN potentials 
show the so-called Coester BandE3 structure (Fig. 2): no prediction is able to re- 
produce the empirical saturation density and energy simultaneously. If the density 
is predicted right the energy is too high, and if the energy is right the density is 
too high. The predictions discussed so far were obtained in the framework of the 
simplest model for the atomic nucleus: point nucleons obeying the nonrelativistic 
Schrodinger equation interact through a static two-body potential that fits low- 
energy NN scattering data and the properties of the deuteron. This conventional 
model for the nucleus is obviously insufficient to quantitatively describe the nuclear 
ground state. It is, therefore, necessary to go beyond the conventional scheme. 

To get some structure into the upcoming considerations, we list in Table 1 
four items relevant to the nuclear many-body problem. The assumptions made 
in the conventional approach are given in column two. Column three states some 
obvious ideas leading beyond the conventional model. Notice that the distinction 
between the various items and assumption is artificial and that there is a great 
deal of overlap among them. Step by step, we will now include most of the possible 
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Table 1. Basic assumptions underlying the nuclear many-body problem 



Simplest assumption Possible extension(s) 



Item 

Degrees of freedom 

Hadron structure 
Interaction(s) 

Dynamical equation 



Nucleons only 

Point structure 

Static, instantaneous 
two-body potential 

Non-relativistic 
Schrodinger equation 



Mesons, 
Isobars, 

Quarks and gluons 

(Quark) sub-structure 

Non-static interactions, 
Many-body forces 

Relativistic 
Dirac equation 



extensions in nuclear matter calculations. One great advantage of Brueckner theory 
is that it allows to take into account additional degrees of freedom in a fairly easy 
and straightforward way. 

2 Meson Degrees of Freedom 

Ever since the meson hypothesis was formulated, it was (at least in principal) 
clear that the full nuclear many-body problem should include nucleons and mesons. 
Nevertheless, traditionally only nucleons have been considered, these interacting via 
a static two-body potential. Even in cases where the two-nucleon force was derived 
from meson theory, the mesons were usually "forgotten" as soon as the nuclear force 
was constructed. "Meson theory" was merely used to provide a suitable ansatz for 
the two-body potential with a convenient parametrization in terms of masses and 
coupling parameters. Thus, the dynamical presence of the mesons was ignored. 
Obviously, from a more fundamental point of view, this is not satisfactory. 

Therefore, we will now include meson degicees of freedom in Brueckner theory 
following the approach suggested by Schiitte.Ej This approach is fieldtheoretic in 
nature, treating baryons and mesons a priori on an equal footing. However, a prin- 
cipal problem of every field-theoretic many-body theory is how to take into account 
the effects of the many-body environment on the particles and their interactions 
(e. g. the single-particle energies in the medium, propagation in the medium, etc.). 
This is difficult to do in a covariant way. Therefore, Schiitte suggested to use 
time-ordered ( "old-fashioned" ) perturbation theoryEj which is similar to the usual 
perturbation theory of ordinary quantum mechanics. Thus, methods familiar from 
non-relativistic many-body theories can be used. 

Starting from a ficld-thcorctic Hamiltonian for mesons and nucleons 

h = t + W, (1) 
where t denotes the unperturbed Hamiltonian (i. e., the operator for the free rela- 
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tivistic energies of nucleons, Em = \J M"^ + q^^^, and mesons, = -^m| + ) and 
W the meson-nucleon interactions, and applying time-ordered perturbation theory, 
the lowest order contribution to free-space two-nucleon scattering is of second order, 
namely. 

Viz) = W W , (2) 

with z the relativistic free energy of the two interacting nucleons (i. e., z = E^+En 
for micleon m and n). Equation (2) ean be understood as an one-boson-exchange 
(OBE) "quasi-potential" (which is energy dependent). The Lippmann-Schwinger 
equation for free-space two-nucleon scattering is 

Tiz) = V{z) + V{z) ^T{z) . (3) 

z — t + le 

Turning to the nuclear many-body problem, we introduce the single-particle 
potential U and rewrite the Hamiltonian Eq. (1) by 

h = ho + hi (4) 

with the unperturbed Hamiltonian 

ho=t + U (5) 

and the perturbation 

hi=W-U. (6) 

As a consequence of this, the one-boson exchange, Eq. (2), is modified in the nuclear 
medium where it reads 

V{z) = W _ \ W. (7) 
z — ho 

Similar to conventional Brueckner theory, the equation for the Brueckner G matrix 
is 

G{-z) = Viz) + V{-z)j^G{z) . (8) 
z — rio 

This equation differs from the (free-space) Lippman-Schwinger eqiiation (3) by the 
Pauli projector Q which projects onto unoccupied two-nucleon states (giving rise to 
the Pauli effect) and by the energy denominator [0 — /lo] which replaces [z — t + ie] of 
Eq. (3) (causing the dispersion effect). Both these effects are the crucial saturation 
mechanisms of conventional Brueckner theory. Note that in all results presented in 
this paper the conventional Brueckner effects are always included. In lowest order, 
the energy per nucleon in nuclear matter is given by 

1 = ^^^™ + ^ E {mn\G{z)\mn-nm)-M (9) 

with M the mass of the free nucleon. The single particle potential used in these 
calculations is defined by 

U{m) =Re {mn\G{z)\mn - nm) (10) 

n<kF 
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Figure 3. Results from nuclear matter Brueckner calculations that include medium effects on 
the two-nucleon interaction. The dotted line is obtained without such medium effects; the dash- 
triple-dot curve includes the medium effect on meson propagation; in addition to this, the dash-dot 
curve contains the dispersion effect from diagrams with intermediate A states. Finally, the dashed 
curve includes all medium effects (i. e., also the Pauli effects from iVA intermediate states). All 
calculations are conducted in the two hole-line approximation using the continuous choice. 



which is appHcd for nucleons below and above the Fermi surface (continuous choice) . 
In the nuclear medium, 

z = e„ + e„ , (11) 

with 

Em = -E'm + , (12) 

in accordance with the unperturbed Hamiltonian ho. 

Formally, Eqs. (8)-(10) are essentially identical to ordinary Brueckner theory. 
The difference is in the quasipotential V which depends on the nuclear medium 
(and its density). The impact of this medium effect is demonstrated in Fig. 3 
where the dash-triple-dot curve is calculated with this effect while the dotted curve 
is without it. It is seen that the medium effect on meson propagation is such that 
the binding energy per nucleon is slightly reduced (in the order of 2 MeV at nuclear 
matter density). This effect is similar to the dispersion effect of ordinary Brueckner 
theory, but, much smaller. The density-dependence of this mesonic effect is such 
that the saturation point moves along the Coester band and not off it. 



3 Isobar Degrees of Freedom 

The excited states of the nucleon play an important role in a genuine and realistic 
meson-exchange model for the nuclear force. The lowest-lying pion-nucleon reso- 
nance, the A(1232) isobar, is essential for NN scattering at low and intermediate 
energies. It provides a large part of the intermediate range attraction and it gener- 
ates most of the inelasticity above pion-production threshold. Again, being aware 
of how crucial this degree of freedom is for the two-body interaction, we should not 
"freeze it out" in the many-body problem. 
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Figure 4. Two- and three-body forces ercated by the A isobar. Sohd hncs represent nuclcoiis, 
double lines A isobars, and dashed lines w and p exchange. Part (a) is a contribution to the 
two-nucleon force. The double slash on the intermediate nuclcon line is to indicate the medium 
modifications (Pauli blocking and dispersion effects) that occur when this diagram is inserted into 
a nuclear many-body environment. Part (b) is a three-nucleon force (3NF). 



The lowest order in which isobars can contribute to the NN interaction is the 
fourth order in the interaction Hamiltonian W (which corresponds to a two-meson 
exchange) . The general structure of the fourth order perturbation is for the case of 
free scattering 

V^'^\z) = W — -W — -W — -W (13) 

^ ^ z-t + ie z-t + ie z-t + ie ^ ' 

with W the meson baryon interaction. The irreducible part of these diagrams 
contributes to the "kernel" V{z) of the scattering equation. One of the diagrams 
generated by Eq. (13) is shown in Fig. 4(a). When inserted into the many-body 
problem the contribution V'^^\z) is altered in a characteristic way; namely it is 
replaced by 

Q Q Q 



yW(j 



w- 



-w- 



-w- 



-w 



(14) 



ho z — ho z — ho 

where the Pauli operator Q projects nucleons onto unoccupied nucleon states. 
One can distiguish between two ways in which the medium exercises influence: 

• The Pauli projector Q cuts out the lower part of the nucleon spectrum in 
intermediate states; this leads to the so-called Pauli effect. 

• The propagator [z — t + ie\~^ is replaced by [z — ho]~^; the effect caused by 
this replacement has become known as dispersion effect. 

Both effects reduce the absolute size of the diagram. Thus, for an attractive diagram 
there is a net repulsive medium effect. Notice that these two effects are similar to 
the conventional Brueckner effects discussed below Eq. (8). 
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Figure 3 includes the results employing the field-theoretic model just sketched. 
The inclusion of the dispersive effects under consideration moves us from the dash- 
triple-dot curve to the dash-dot curve, and the Pauli efi^ects bring us further up to 
the dashed curve. This group of medium effects are substantially larger than those 
from the OBE terms discussed in the previous section. 

The results niesented in Fig. 3 are based upon the Bonn Full Model for the 
NN interaction.cEl This model includes also the (repulsive) np diagrams (for which 
the medium effect causes a net attraction). However, since the sum of 2t: and Trp 
diagrams is attractive, the 27r-exchange being dominant at intermediate range, the 
net medium effect is repulsive. It is clearly seen that dispersion and Pauli effects are 
about equally important, the latter typically increasing more strongly with density. 
We note that these calculations also include all non-iterative diagrams (stretched 
and crossed box diagrams) involving A isobars which contribute about as much to 
the medium effects as the box (iterative) diagrams. 

The density-dependence of the effects due to A degrees of freedom is only slightly 
stronger than that of the conventional saturation mechanisms, bringing the satu- 
ration point not markedly off the Coester band. 



4 Many-Body Forces 

When discussing many-body forces, caution is in place. Notice that the diagram 
Fig. 4(b) is a genuine three-body force in a model space which consists of nucleon 
states only. In an extended Hilbert space which includes A isobar states. Fig. 4(b) 
represents just a three-particle correlation. Thus, most many-body forces are arti- 
ficially created by freezing out degrees of freedom. In this respect they are merely 
artefacts of the particular theoretical framework applied. Because of this model 
dependence of the terminology, it is useful to introduce an operating definition for 
n- nucleon forces which we will take to be the following: forces that depend in an 
irreducible way on the coordinates or momenta of n nucleons when only nucleon 
degrees of freedom are taken into account. 

The fact that most many-body forces are an artifact of theory implies the follow- 
ing rule for how to deal with the many-body-force issue in a proper way: when you 
introduce a new degree of freedom, take it into account in the two- and many-body 
problem, consistently. Field-theoretic models for the 27r-exchange contribution to 
the NN interaction require the A(1232) isobar, which also creates n-nucleon forces 
in the many-body system. Consequently, when introducing the isobar degree of 
freedom, it should be included in the two- and many-nucleon forces simultaneously. 
In Fig. 4(a) we showed a contribution to the two-nucleon interaction that involves 
a A isobar in the intermediate state; the existence of the A isobar implies that, in 
the many-body system, the diagram Fig. 4(b) will occur (and many other diagrams 
involving As) which — according to the above definition — represents a three-nucleon 
force (3NF). 

For a spin- or isospin saturated system, the lowest order term. Fig. 4(b), van- 
ishes. However, when a third interaction between the nucleon legs is introduced, 
then there is a significant attractive contribution. A very thorough calculation p£ 
diagrams of this kind has been conducted by Dickhoff, Faessler, and Miither.Ell 



mbx: submitted to World Scientific on February 9, 2008 





^ -10 
> 

<; -20 
-30 



1 1.2 1.4 1.6 

kp (fm-') 

Figure 5. Many-nucleon force contributions to the nuclear matter energy. The dashed Hne is 
identical to the dashed line of Fig. 3. Adding the 3rd and 4th order ring diagrams with A- 
isobars results in the dash-dot curve. The latter contributions are from the Tiibingen group'^^ 
who produced results for kp < 1.6 fm~^ only. Further contributions from meson and nucleon 
renormalization lead to the solid curve which is our final result. The dotted curve is identical 
to the dotted curve of Fig. 3 and includes neither medium effects on the two-nucleon force nor 
many-nucleon force contributions. 



Speaking in terms of the hole-line expansion, these authors have caculated the 
three and four hole-hne contributions of the ring type with isobar degrees of free- 
dom. In terms of the 'nucleons only' language, the ring diagrams involving isobars 
are ccpiributions from many-nucleon forces. The results obtained by the Tiibingen 
groupEZI are shown in Fig. 5. 

The top line (dashed) in Fig. 5, which is the starting point of our quantitative 
considerations in this section, repeats the final result of the previous section; it 
is obtained in lowest order Brueckner theory (two hole-line approximation) and 
includes all medium effects on the two-nucleon interaction discussed in Sects. 2 
and 3. Now, we add the attractive contributions from ring diagrams of third and 
fourth order involving A isobars (i. e., 3NF and 4NF contributions) and arrive at 
the dash-dot curve in Fig. 5. 

Another class of diagrams that belongs into the category jaf many-nucleon forces 
(according to the above definition) was calculated in Ref. E3 where it was shown 
that the selfenergy diagrams of tt and p together with nucleon selfenergy corrections 
result in a net attractive effect (which in terms of magnitude is of the same size as 
the repulsive mesonic effect discussed in Sect. 2). When we add this contribution 
to the dash-dot curve in Fig. 5, then we obtain the solid curve which is our final 
result. In this figure, we also show again the dotted curve of Fig. 3 which represents 
the result when no medium effects and no many-nucleon forces were included. The 
comparison of the dotted curve with the solid curve in Fig. 5 reveals that the sum 
of all medium effects on the NN force plus all many-nucleon force contributions 
is much smaller than the individual contributions; or in other words, there are 
large cancelations. However, the cancelations are not perfect and the net effect is 
moderately repulsive. 

To summarize this and the previous section, isobar degrees of freedom have 
essentially two consequences in nuclear matter: 
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• medium effects on the two-nuclcon interaction and 

• many-nucleon force contributions. 

Both are large effects/contributions — but, of opposite sign. In a consistent treat- 
ment of degree(s) of freedom either both effects occur simultaneously, or none. One 
of these two effects alone, in isolation, does not exist in reality. Therefore, to take 
into account only one of them (for instance, only the three-body force contribu- 
tions, ignoring the medium effects on the corresponding two-body diagrams) yields 
a substantially distorted picture. In fact, the almost cancelation between these two 
effects/contributions may be the deeper reason why, ultimately, many-body forces 
may not play a great role in nuclear physics; it may also be the reason why the 
traditional two-body force picture is by and large rather successful. 

We note that our findigs in nuclear matter are supported by accurate and inde- 
pendent calculations conducted for the triton binding energy. Triton calculations 
of this kind were first performed by the Hannover group.tJ Recently, these calcula- 
tions have been improjzBd and extended by Picklesimer and coworkers^ using the 
Argonne V28 A-modelB3 They find an attractive contribution to the triton energy 
of -0.66 MeV from the 3NF diagrams created by the A, and a repulsive contribu- 
tion of -1-1.08 MeV from the dispersive effect on tlie two-nucleon force involving A 
isobars. The total result is 0.4^2 MeV repulsion^ This means that — also in the 
triton — the final result is a moderate repulsion, similar to our findings in nuclear 
matter. The consistency of the results for the two very different many-nucleon sys- 
tems suggests that there may be some general validity and model-independence to 
our findings. 

Besides the A-isobar, there are other mechnisms that can give rise to three- 
nucleon forces. In recent years, it has become fashionable to consider chiral La- 
grangians for nuclear interactions. Such Langrangians may create diagrams which 
represent effective thre&-nucleon potentials. However, Weinberg has shown that all 
these diagrams cancelfiJ Moreover, based upon power countiag arguments, one can 
show that chiral n-nucleon forces with n > 4 are negligible.c3 

5 Comparison to Simple Models 

In Sections 2 to 4, we have presented the results of the most comprehensive and con- 
sistent nuclear matter Brueckner calculations ever conducted in terms of hadronic 
degrees of freedom — with the final result given by the solid curve in Fig. 5. After all 
this work which took the research groups involved about a decade, one may raise a 
very naive, but reasonable question: how does this final result compare with the pre- 
dictions from the usual simplistic approach in which a static (energy-independent) 
two-nucleon potential is applied? 

To answer this question, we show in Fig. 6 the saturation curves generated by 
some NN potentials together with the 'final result' of Sect. 4 (thick solid curve, 
identical to solid curve of Fig. 5; note also that Figs. 5 and 6 have slightly dif- 
ferent scales). The NN potentials applied in Fig. 6 are all fropLlha-new family 
of high-precision/high-quality potentials developed in the 1990's.OE3c3 These new 
potentials fit the NN data below 350 MeV laboratory energy with the 'perfect' 
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Figure 6. Nuclear matter predictions by some recent two-nucleon potentials. The dashed, dash- 
dot, and dash-triple-dot curves are produced by the (local) Argonne Vig,*'* Nijm-II,*^ and 
ReidgS**^ potentials, respectively; and the thin solid line is predicted by the (nonlocal) CD-Bonn"*^ 
potential. The thick solid curve is identical to the solid curve of Fig. 5 and is shown for com- 
parison. Note that this curve is given only for kp < 1.6 fm"'^, where it cannot be distinguished 
from the CD-Bonn curve on the scale of the figure. The star is the saturation point predicted by 
Brueckner and Gammel^ in 1958. All calculations are based upon the two hole-line approximation 
using the continuous choice. 



/datum « 1; consequently, all these new potentials produce essentially identical 
predictions for the two-nucleon observables. However, this does not imply that these 
potentials are identical. Indeed, they differ in their off-shell behavior as indicated 
by the fact that some_C|f these potentials are local and some nonlocal. 

The Areonne Vis^ Nijm-II,El and Reid93El potentials are local and their nu- 
clear matter predictions are represented in Fig. 6 by the dashed, dash-dot, and 
dash-triple-dot curves, respectively. The CD-Bonn potential is nonlocal and pre- 
dicts the thin solid curve. 

Figure 6 clearly reveals that all local potentials lead to almost identical nuclear 
matter results. This statement is even true if one includes the result by Brueckner 
and Gammel of 1958 with the local Gammel-Thaler potential (see star in Fig. 6). 
So, in terms of local potentials, there has not been much progress in the past 40 
years. 

The story is different with the nonlocal CD-Bonn potential which predicts a 
distinctly different result. Moreover, the nonlocal CD-Bonn reproduces amazingly 
wellc3 the final result from the thorough and comprehensive calculations of Sects. 2 
to 4. 

Thus, the bottom line is that the result of the 'faithful' calculation that takes 
meson, nucleon and A degrees of freedom consistently into account can, indeed, be 
reproduced by a two-nucleon potential. However, not by any potential, it has to be 
a nonlocal one. Or in other words, the nonlocality, which is one important element 
of the comprehensive theory of Sects. 2 to 4, cannot be ignored. 

Since the CD-Bonn is based upon relativistic meson exchange, it contains the 
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same nonlocalities as the fieldtheoretic models applied in Sects. 2 to 4, except for the 
energy-dependence. This energy-dependence of the NN interaction can obviously 
be neglected — the reason being the large cancelations between the medium effects 
on the two-nucleon interaction (which are caused by the energy-dependence) and 
the n-nucleon force contributions. However, the other forms of nonlocality cannot 
be neglected. In short: you can make it simple but not too simple. 



6 Relativistic 'Dirac' Effects 



We will now explore an important effect not included in the systematic investiga- 
tions conducted up to Sect. 4. Notice that, in the fieldtheoretic approach sketched 
in Sect. 2, one can go one step further and use the Dirac equation for single-particle 
motion in nuclear matter 

(^ - A/ - f7)u(p, s) = (15) 

or in Hamiltonian form 

{a • p + PM + pU)u{p,s) =^ epu{p,s) (16) 

with 

U = Us + fUv (17) 

where Us is axL attractive scalar and Uv the time component of a repulsive vector 
field. (Notation as in Ref. (5 = 7°, a' = etc.) The fields, Us and Uv, are 

in the order of several hundred MeV and strongly density dependent. The solution 
of Eq. (15) is 



with 




"(P' = \ ',~r ^LlP. Xs (18) 



M^M + Us, Ep = ]VP + p-^ , (19) 

and Xs a Pauli spinor. 

The Brueckner equation now reads 

G{z) = V + V^^G{S). (20) 
z - ho 

The essential difference to standard Brueckner theory is the use of the potential V. 
As indicated by the tilde, this meson-theoretic potential is evaluated by using the 
spinors Eq. (18) instead of the free Dirac spinors applied in scattering (and con- 
ventional Brueckner theory). Since Us (and M) are strongly density dependent, so 
is the potential V. M decreases with density. The essential effect in nuclear mat- 
ter is a suppression of the (attractive) a-exchange; this suppression increases with 
density, providing additional saturation. It turns out (see Fig. 7 below) that this 
effect is so strongly density-dependent that the empirical nuclear matter saturation 
energy and density can be reproduced simultaneously. 
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The lowest-order energy of nuclear matter is 

(21) 

The single-particle potential, 

U{m) = ^(TO|[/|m) = ^{m\Us + 7°t/y|m) = ^C/g + Uy , (22) 
is determined from the G-niatrix in formally the usual way, 

U {in) — Re — ; — —{mn\G{z)\mn ~ mn) (23) 

n<fcF EnE„i 

with m below and above the Fermi surface (continuous choice) and 

2 = e-m + En , (24) 

where the single particle energy is given by 

M 

em= -^{m\'i •Pr„, + M\m) + U{m) (25) 

Em 

^ Era + Uv ■ (26) 

Note that, in this approach, the nucleon states |m) and \n) are represented by 
Dirac spinors of the kind Eq. (18) and an appropriate isospin wavefunction, {m\ and 



(n| are the adjoint Dirac spinors u — u^^^; Em = y -^^^ Pm- The normalization 
of the Dirac spinors is um = 1. 

Results are shown in Fig. 7 where the lower solid and lower dashed curves repeat 
the predictions presented in the previous section for the Argonne Vls and the CD- 
Bonn potentials, respectively. To these curves, we now add the relativislir; Dirac 
effects based upon the formalism just sketched and calculated in Ref. c3. This 
results in the upper dashed and solid curves for Vig and CD-Bonn, respectively. 

As discussed, this Dirac effect allows a correct reproduction of the saturation 
density and energy — as demonstrated by the upper CD-Bonn curve. It is also 
evident that slight overbinding in the nonrelativistic approach is a necessary pre- 
requisite for ultimately predicting density and energy right. As revealed by the 
upper dashed line in Fig. 7, local potentials lack binding energy and, therefore, 
cannot reproduce nuclear saturation correctly. 

Alternatively, one may take a phenomenological approach in which one tries 
to simulate the. Dirac effect by a repulsive 3NF. This method is pursued by the 
Urbana groupll3 and we have included their results in Fig. 7 in terms of the dotted 
lines. The lower dotted curve is the variational result with the Vig two-nucleon 
potential and should be compared with the lower dashed curve (Brueckner with 
Vis). The closeness of these two curves confirms what we discussed in conjunction 
with Fig. 1, namely, a satisfactory agreement between Brueckner theory and the 
variational approach. The Urbana group developed a repulsive phenomenological 
3NF designed to improve the nuclear saturation prediction. Adding this 3NF leads 
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Figure 7. The impact of relativistic Dirac effects on nuclear matter saturation. The lower solid 
and dashed curves are without and the corresponding upper curves are with Dirac effects. Dashed 
curves represent Brueckner results based upon the Argonne Vig potential and solid curves corre- 
spond to CD-Bonn. Dotted curves are variational results with Vig where the upper dotted curve 
includes a repulsive 3NF (from Akmal et al.'^^). 

US from the lower dotted curve to the upper dotted curve. The latter is obviously 
close to the upper dashed line which implies that the repulsive 3NF is simulating 
the Dirac effects well. The upper dotted and dashed curves are the final results for 
Argonne Vis and prove that — no matter what mechanisms are invoked to get the 
saturation density right — local potentials always underbind nuclear matter. 

In conclusion, nonlocality of the nuclear force is crucial to explain nuclear bind- 
ing quantitatively. This should be no surprise since the nuclear force is inherently 
nonlocal, anyhow. 

7 Conclusions 

Using Brueckner theory, we have investigated systematically the influence of pfipn- 
nucleonic degrees of freedom on the predicted properties of nuclear matter .e2I In 
general, many-body calculations of this kind can be very involved and most cum- 
bersome. However, in the framework of Brueckner theory, these rather sophisticated 
investigations can be conducted in a straightforward and managable way. This is 
one of the great advantages of Brueckner theory. 

One result of our investigation is that explicit meson and A-isobar degrees of 
freedom do not improve the saturation properties of nuclear matter. However, they 
teach us something about many-body forces: a mechnism that has the potential of 
creating many-nucleon forces should be treated consistently in the two- and many- 
nucleon system. Even though individual contributions caused by this mechanism 
can be large, the net result turns out to be very small — due to characteristic cance- 
lations between effects from two- and many-nucleon terms. Indeed, the final result 
is very similar to what is obtained from a nonlocal two-nucleon force, since nonlocal- 
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ity is an inherent quality of the more comprehensive and fundamental fieldtheoretic 
approach. If the force is local, insufficient binding is predicted. 

Wc also discussed an extension of Brucckncr theory in which the rclativistic 
Dirac equation is used for single particle motion (also known as Dirac-Brueckner- 
Hartree-Fock). This has a substantial impact on nuclear saturation such that the 
empirical density and energy can be reproduced. 

Based upon the work summarized in this talk and the countless papers not 
mentioned, we conclude that Brueckner theory has contributed in a crucial way to 
our understanding of the nuclear many-body problem. 
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